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Abstract. We apply tomography, a general method for reconstructing 3-D 
distributions from multiple projections, to reconstruct the momentum distribution 
of electrons produced via strong field photoionization. The projections are 
obtained by rotating the electron distribution via the polarization of the ionizing 
laser beam and recording a momentum spectrum at each angle with a 2-D 
velocity map imaging spectrometer. For linearly polarized light the tomographic 
reconstruction agrees with the distribution obtained using an Abel inversion. 
Electron tomography, which can be applied to any polarization, will simplify the 
technology of electron imaging. The method can be directly generalized to other 
charged particles. 



PACS numbers: 32.80.Rm, 33.60.+q 

1. Introduction 

Photoelectron spectroscopy is an important method to study atoms [1 , molecules 
[2] and condensed matter [3]. The photoelectron energy, momentum and angular 
distributions can identify the state of the residual ion provide insight into 
correlation effects [5], or yield structural information [6l [7]- Time resolved 
photoelectron spectroscopy allows us to study the dynamics of these observables in 
evolving systems [H]. The technology for measuring electrons is very well advanced. 
Single electrons can be detected with ^50% efficiency using electron multiplier 
tubes and microchannel plates (MCPs). To measure the three dimensional electron 
momentum distribution a number of approaches have been developed. For example, 
an array of single channel electron multipliers can obtain two dimensional information. 
The time-of- flight then adds the third dimension [9]. A single detector can replace 
the array if the electron is captured in an electric and magnetic field and recorded 
with a combination of MCP and position sensitive detector. When used with a 
multihit capable delayline detector, the complete 3-D momentum distribution of single 
electrons can be measured |10| . However, this powerful method is experimentally 
demanding and comes at the cost of reduced resolution. Furthermore, the 3-D 
momentum of every single electron is not always necessary since the accumulated 
electron distribution in the laboratory frame is often a sufficient observable. 
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If an MCP is used in conjunction with a phosphor screen and a CCD camera, the 
technological side is greatly simplified, and provides, in continuous operation, easy 
access to the lab frame distributions of electrons and ions. In a widely employed 
detection scheme known as velocity map imaging (VMI) [TT], only electrostatic fields 
are used to guide the electrons to the detector. Through an integrated electrostatic lens 
a direct mapping of the initial photoelectron velocity distribution onto the detector is 
achieved. To obtain a An collection angle these electric fields are typically several kV in 
strength resulting in time-of-flight distributions on the order of the temporal response 
of the detector. Therefore, in VMI the information along the axis perpendicular to 
the detector is lost. 

Previous approaches to determine the 3-D momentum distribution in VMI 
experiments used inverse Abel transforms [121 1131 E], slice imaging [T^, or time- 
resolved event counting [HI [17]. Each of these techniques suffers limitations when 
applied to non-symmetric distributions. By construction the Abel methods are limited 
to cases with cylindrical symmetry. The other methods can be used to image non- 
symmetric distributions in 3-D provided the particles are suitably spread in time-of- 
flight. For electrons this is difficult and new approaches are needed. 

A recent attempt to address this problem was to use a tomographic method to 
reconstruct the momentum distribution in three dimensions |181 119j . This allows the 
reconstruction of a non-symmetric distribution in a high data throughput experiment 
like VMI. 

We build on the idea of WoUcnhaupt et al [18] and show how this technique can 
be applied to tunnel ionization of an atom. In experiments with strong laser fields, 
control of the polarization is essential for affecting electron re-collision processes on 
attosecond time scales. Imaging the 3-D photoelectron distribution in arbitrarily 
polarized light therefore provides useful information for strong field science. We 
apply the filtered backprojection technique widely used in medical CT to reconstruct 
the electron distribution. This is in contrast to WoUenhaupt et al who used a 
Fourier approach. The filtered backprojection technique can give exact reconstructions 
compared to approximate results obtained with unfiltered techniques |20j . Using 
multiphoton ionization of argon as an example, wc make a quantitative comparison 
between the tomographic method and the Abel inversion method for the case where 
an Abel inversion is possible. We present the reconstructed 3-D electron momentum 
distributions created in hnearly and elliptically {E1/E2 = 0.89) polarized light in the 
tunnelling regime. 

2. Methods 

2.1. Velocity map imaging 

In velocity map imaging experiments an inhomogeneous electric field is used to 
project the velocity vector of charged particles into the 2-D plane of the detector. 
This is shown schematically in figure [Ta| The electric field acts as a lens which 
'pancakes' the 3-D velocity distribution, preserving it's magnitude along the x' axis 
and integrating over the y' axis. The vectors vi,V2,V3 all having the same Vx' but 
different Vy' are projected onto the same spot on the detector. For this to work it 
is necessary that the initial velocity of the particles is negligible compared to their 
acceleration in the spectrometer's DC field. The y' component is not observable 
unless additional information is provided, such as symmetry assumptions or, as we 




(a) Principle of velocity map imaging. (b) Principle of tomographic imaging. 
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(c) Experimental setup for VMI tomography. 



Figure 1: (a) The principle of VMI. An electrostatic lens projects the x' component 
of the velocity vector (dark arrows) into the plane of the detector, (b) The principle 
of tomographic imaging. A distribution / {x, y) is shown at two different orientations 
and projected onto the detector along the electric field E. (c) Experimental setup for 
VMI tomography. See text for details. 

will show, measurements from multiple directions. The former is a condition for Abel 
inversion and allows the inversion of a single projection. The latter is a requirement 
for tomographic inversion and is a more general approach. 

2.2. Tomographic imaging 

The principle behind tomographic imaging [21j is sketched in figure [Tbl An arbitrary 
distribution / (x, y) is shown at two different orientations relative to a detector (shown 
as the horizontal lines). The detector frame of reference {x',y',z') is related to the 
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distribution frame of reference by, 
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For the application of the tomographic method the rotated distribution must be 
projected onto the detector. Mathematically this is given by the Radon transform, 



Physically, the Radon transform is carried out in a VMI experiment by the 
spectrometer's DC electric field, oriented along the long arrows in figure llbl The 
projection of the velocity vector discussed in section 12.11 is the experimental Radon 
transform. Examples of the projections of / (x, y) at two different orientations are 
sketched as the curves Pq {x') and Pgo {x') in figure (Tb] By collecting projections 
for many different orientations the tomographic algorithm is able to reconstruct 
the complete velocity distribution. Unlike medical tomography, wherein a detector 
typically covers a portion of the volume of interest, the projections in VMI experiments 
span the complete momentum space of the photoelectrons. Thus, by rotating the 
distribution there is enough information to recover the full 3-D distribution for any 
polarization state. 

2.3. VMI tomography 

The setup for VMI tomography is depicted in figure [Tel The laser propagates along 
the z' {z) axis with polarization in the x'y' (xy) plane. In figure [Tel the polarization is 
shown as linear, but in general it can be elliptical. A half wave plate is used to rotate 
the axes of the polarization. The ability to rotate the polarization and thereby the 
velocity distribution in the chamber is a necessary step in tomographic reconstruction. 
In section [3l results are presented for linear and elliptical polarization. 

The beam then enters the VMI spectrometer where it is focused to the centre 
of the spectrometer using a parabolic mirror (not shown in figure ITcl) . The laser 
focus is orthogonally crossed by a supersonic gas jet having a density corresponding 
to 5 X 10~^ torr. The VMI electrodes are negatively biased with the correct voltage 
ratio to project the electron velocity distribution (shown as the dark cloud in figure 
[Tell onto the 2-D detector. The electrons are accelerated to roughly 2.5 keV in the 
spectrometer's DC field or about 140 times their kinetic energy acquired from the 
laser pulse. The detector is pulsed by a high voltage switch to eliminate noise and 
dark signal. The timing gate on the high voltage pulse is 300 ns long. The amplified 
signal from the detector is read by a CCD camera acquiring continuously at 12 fps. 
The Ti:Sa laser system used for these experiments produces 50 fs pulses at 500 Hz 
centred at 805 nm with a bandwidth of 22 nm. 

2.4- Image inversion 

The inversion of the 2-D projections follows the parallel ray filtered backprojection 
algorithm [21j . This method can yield exact reconstructions in contrast to techniques 
without a filter step. In any computer implementation of the algorithm, the projections 
are always sampled at some distance interval r and some angular interval A9. The 
first step in the tomographic algorithm is to filter each measured projection Pg (x', z). 




(2) 
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This is implemented as the convolution of Pg {x') with a filter function. We have used 
the discrete "Shepp - Logan" filter pO] : 

M-') = - (4;^) (3) 

where x' = nr with r as the sampling interval and n is an integer. The filtered 
projection is 

Qe{x',z) = hix')^Pgix',z) (4) 

As explained in section 12.11 the data were originally acquired by integration along 
the spectrometer's electric field {y' axis). At this stage in the reconstruction, each 
filtered projection Qe (x' , z) must be backprojected along the y' axis. Summing over 
all the projection angles yields the reconstructed distribution. Mathematically, the 
reconstructed distribution F (x, y, z) is given by 

F{x,y,z) - / Qe{x\z) Ad (5) 

K 

= j^'^Qsi {x cos 9i-y sin 9i,z) (6) 
1=1 

where the projections are sampled at K different angles and < 0^ < tt. Implementing 
the discrete backprojection ([6]) requires one-dimensional interpolation to find the 
projection value Qe^ {x') at each point {x^y) in the distribution frame of reference. 

Equation ([6]) is the reconstructed momentum distribution in three dimensions. 
Thus, by rotating the laser polarization the tomographic method allows us to 
reconstruct the complete momentum distribution. 

3. Results 

First, we compare the tomographic method with the inverse Abel method for linearly 
polarized light. For an Abel inversion the polarization axis must be parallel to the 
detector face. A single two-dimensional projection was recorded at a laser intensity 
of 2 X 10^^ W/cm^ and then inverted using the inverse Abel method. A profile along 
the polarization axis in the inverted distribution is shown as the solid line in figure 
[2] The profile has been normalized to its maximum value. The distribution shows 
multiple peaks spaced by one photon's energy. This is characteristic of above threshold 
ionization [22] . 

Next, the tomographic method was tested under the same conditions. The laser 
parameters were held constant and the polarization direction was rotated using a 
broadband half-waveplate (Bernhard Halle GmbH). Examples of two projections taken 
at and 90 degrees are shown in figure [2^ and ^p. In figure [2^ the polarization is 
parallel to the detector face, and in figure [lb it is perpendicular. The two dimensional 
projections Pg {x' , z) were recorded as the polarization rotated from — 90 degrees in 
steps of 2 degrees. These projections were then duplicated to yield the tt projection 
angles required by the parallel ray backprojection algorithm. 

The tomographic method was used to reconstruct the momentum distribution in 
15 different slices through the 3-D momentum distribution. A spacing of 0.071 atomic 
units was used between each adjacent slice along the laser propagation (z) axis. The 
result is the 3-D electron momentum distribution. A profile along the polarization 
axis and in the plane of polarization [z = 0) is shown as the dashed line in figure [2] 
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Figure 2: Profile along the polarization axis of the 3-D momentum distribution. Blue 
(solid) line: distribution inverted using the inverse Abel (BASEX) method 14J. Red 
(dashed) line: distribution inverted using the tomographic method. Inset: examples 
of projections used in the tomographic reconstruction a) 6 — 0° , h) 9 — 90° 



Again, the profile has been normalized to its maximum value. The profile generated by 
the tomographic method is reconstructed using a series of 1-D projections measured 
at different angles. This is in contrast to the inverse Abel method which uses a 
2-D projection to infer the structure of the momentum distribution. There is good 
agreement between the tomographic method and the result from the BASEX method 
[14] . This shows the tomographic method reproduces the result of the inverse Abel 
transform for linearly polarized light. 

The power of the tomographic method lies in its application to distributions 
that need not contain any symmetry. A quarter wave plate was inserted into the 
beam and rotated by 43° to create elliptically polarized light at an intensity of 
3 X 10"'^'' W/cm^. Using the half wave plate in figure [Tel the polarization ellipse was 
rotated from — 178 degrees in steps of 2 degrees. At each angle a projection of the 
electron velocity distribution containing roughly 6 x 10^ electrons was acquired. The 
tomographic method was used as before to reconstruct the 3-D velocity distribution. 
The distribution in the plane of polarization is shown in figure [3al 

Ellipticity appears in two ways in the reconstructed image. The most obvious 
evidence is in the nonuniform momentum distribution. This is shown more clearly in 
figure [3b] where a profile along the major axis is shown as the red line; the minor axis 
is shown as the black line. In figure [3b] it is clear that there are more electrons along 
the minor axis of the ellipse. This counterintuitive result is explained in reference 
[23] , It occurs because the ionization rate is a function of the laser's instantaneous 
electric field while the drift velocity is a function of the vector potential at the phase 
of ionization. Since the two are 7r/2 out of phase, the larger ionization rate appears 
along the minor axis of the velocity distribution, as shown in figure [3b] Using the size 
of each of the axes, the laser ellipticity is found to be 0.89. 
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(a) (b) 



Figure 3: Tomographic imaging of the electron velocity spectrum in elliptically 
polarized light, (a) Reconstructed distribution in the plane of polarization. The 
colour code represents the relative number of electrons, (b) The number of electrons 
as a function of their speed measured along the minor (black) and major (red) axes. 




(a) linear polarization (b) elliptical polarization 



Figure 4: 3-D reconstructions using the tomographic method. The colour code is the 
same as figure |3al Values less than 0.1 are transparent. 

The tomographic method can be applied to visualize the 3-D velocity distribution 
for any polarization state. Figure 3] shows a three-dimensional representation of the 
reconstructed velocity distribution in both linearly and elliptically polarized light. 
A surface of constant probability is shown in navy blue. On the interior of the 
distributions the colour coding reflects the probability to measure an electron at each 
point in the 3-D distribution. In figure [4] the data was averaged over neighbouring 
volume elements. 
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4. Conclusion 

We have shown how ideas of tomographic imaging can be extended to study tunnel 
ionization in velocity map imaging experiments. No assumptions are needed about 
the symmetry of the measured distribution. The method accurately reproduces the 
3-D momentum distribution in linearly polarized light. We have applied the method 
to reconstruct the velocity distribution of electrons in multiphoton ionization of argon 
by elliptically polarized light ~ a case not handled by inversion techniques based on 
the Abel transform. 

There have been some recent papers [231 [3S] questioning the applicability of 
tunnelling models 26J to strong field phenomena in 800 nm or even 1400 nm laser light. 
Experiments with circular or elliptical polarization will allow the tunnelling model to 
be tested without the complexity of re-collision. Resolving the ionized electrons in 
three dimensions will give new insight into the ionization mechanism. 
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